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Q Abstract. The models developed to describe the spectral energy distribution (SED) of blazars 

can be divided into leptonic or hadronic scenarios, according to the particles responsible for the 
high-energy component. We have developed a new stationary code which computes all the relevant 
leptonic and hadronic processes, permitting the study of both leptonic and hadronic scenarios in a 
consistent way. Interestingly, mixed lepto-hadronic scenarios (in which both components contribute 

^ to the high energy emission) naturally arise in this framework. We present the first application to 

pq the well known BL Lac object PKS 2155-304. 
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2 INTRODUCTION 

, In the framework of the unified active galactic nuclei (AGN) model, blazars are con- 

sidered as radio-loud AGN whose relativistic jet is pointed close to the line-of-sight 
[1]. Their SED is characterised by two bumps, one peaking at low energies (infrared to 
^ X-rays), the other in y-rays. While the first component is unanimously ascribed to syn- 

<N chrotron emission by a non-thermal population of electrons and positrons, the origin of 

^ the second component is still under discussion. In leptonic models it is ascribed to an 

inverse Compton process between the electrons and a low energy photon field (their own 
2 synchrotron emission, or external photons), while in hadronic models it originates from 

CS| synchrotron emission by protons and secondary particles coming from p-y interactions. 

We have developed a new stationary lepto-hadronic code which describes all the rele- 
> vant process, permitting the study of both scenarios in a consistent framework, as well 

^ as interesting mixed lepto-hadronic scenarios. We first provide the details of the code, 

;-i and we then present an application to the BL Lac object PKS 2155-304. 



DESCRIPTION OF THE CODE 



The code we present is a direct evolution of the code developped by Katarzyhski et al. 
[2] Several modifications and improvements have been applied to the original leptonic 
code: more details can be found in [3]. 

The emitting region is assumed to be filled with a stationary population of relativistic 
electrons/positrons and protons. The proton energy distribution Np{yp) is described by a 



power-law function with index Up, defined above Jp-^min and with an exponential cut-off 

at Yp;Max- 

The proton population in the emitting region interacts with the low energy photons 
through the photo-meson process and the electron-positron pair production (Bethe- 
Heitler process). The photo-meson interaction is evaluated using the publicly available 
Monte-Carlo code SOPHIA [4] which gives as output the distributions of the stable 
particles produced in the interaction. The low-energy photon field is represented by the 
synchrotron emission from the primary electron population in the jet, which represents 
by far the dominant component at low energies for high-frequency-peaked BL Lac 
objects as PKS 2155-304. The most important features of the computation of hadronic 
processes are the following: 

• We first compute the synchrotron emission by protons. 

• Before calling SOPHIA, the energy of the proton is modified in order to take into 
account the synchrotron losses following Miicke and Protheroe [5]. 

• Photons coming from the decay and the synchrotron emission from coming 
from the K"^ channel can trigger an electro-magnetic cascade, supported by syn- 
chrotron emission. The stationary state of the cascade emission has been evaluated 
iteratively computing for each generation the stationary distribution of pairs and 
the associated synchrotron emission. 

• Within SOPHIA, we corrected for synchrotron losses the energy of kaons, pions 
and muons before decay. 

• Within SOPHIA, we retrieved the muon (jl^) spectra before their decay into elec- 
trons and positrons. The stationary spectrum is computed taking into account both 
decay and synchrotron cooling, and their synchrotron emission is then evaluated. 

We consider that the proton index is equal to the electron index before the break (under 
the assumption that the acceleration mechanism is similar), and we fix the maximum 
proton energy through a comparison between acceleration and cooling time-scale (see 
Miicke and Protheroe [5]). With respect to a pure leptonic model, we add then only one 
free parameter in the code, the proton normalization factor. 

MODELLING OF PKS 2155-304 

As a first application we study the BL Lac object PKS 2155-304, using in particular the 
data from the multi wavelength campaign of 2008 [6], in which the source was found to 
be in a low flux state. 

We model the source in four different scenarios, shown in Fig. 1 (the corresponding 
model parameters are given in Table 1): a standard synchrotron-self-Compton model, 
a hadronic model and two different lepto-hadronic models. The Doppler factor of the 
emitting region has been fixed at 30 for all the models. For the SSC model we adopted 
a low magnetic field (B = 65 mG), and a relatively large emitting region (R ~ lO^^cm). 
In the hadronic model, in order to have a significant synchrotron emission from the 
protons, we need to increase the magnetic field up to 80 G, reducing the emitting 
region by two order of magnitudes (R 10^^ cm) in order to fit the low energy bump 



FIGURE 1. Top left: modelling of PKS 2155-304 in a synchrotron-self-Compton scenario. Top right: 
modelling of PKS 2155-304 in a hadronic scenario. Bottom: modelling of PKS 2155-304 in a lepto- 
hadronic scenario. The colour-code is as follow: red-primary electron synchrotron emission; green-inverse 
Compton emission; blue-proton synchrotron emission; light blue-muon synchrotron emission; yellow- 
photons from decay; pink-synchrotron emission from secondary coming from tt* decay. The 
emission from the associated electro-magnetic cascade is shown by the dotted curves (same colour). 

with the primary electron synchrotron emission. Interestingly, the synchrotron emission 
from secondary particles (ji and cascades associated with 7 and e^) coming from p-j 
interactions significantly contributes to the high-energy bump. 

The first lepto-hadronic model uses the same parameters that were adopted for the SSC 
model, but we add protons in the emitting region. Given the low magnetic field, they 
do not contribute directly to the SED, but the emission from secondary pairs can play 
a role in the TeV energy range (in this case the power of the emitting region increases, 
reaching the Eddington luminosity of the source). The second lepto-hadronic model lies 
in between the SSC and the hadronic one, with a magnetic field of 5 G. In this case the 
SSC component dominates only up to a few 10 GeV energies, while at higher energies 
the synchrotron emission by secondary pairs becomes the dominant process. In this case 
the Fermi and H.E.S.S. data would then be ascribed to two different processes. 
An interesting result is that the presence of hadrons affects the y-ray emission, inducing 
a hardening of the TeV spectrum (that can be seen in both the hadronic model and the 
first lepto-hadronic model). This effect can be tested by the next-generation Cherenkov 
telescope CTA. Another important feature is that the cascade component tends to fill 
the gap between the two bumps much more than in pure SSC models. The upcoming 
hard-X-ray satellites NuSTAR and Astro-H will put new constraints on this part of the 
spectrum. 



TABLE 1. Parameters used for the SSC and the hadronic modelling 
of PKS 2155-304. The luminosity of the source has been calculated as 
Ljgt = nR^i.gCr^^ii,{uB + Mg + Up). For the SSC scenario, we included the 
term Up^cold^ which represents the contribution of cold protons in the 
emitting region, estimated following Sikora et al. [7]. 
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